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TECHNICAL MEMORANDUM 


THE EFFECT OF INDUCED CHARGES ON LOW-ENERGY PARTICLE 
TRAJECTORIES NEAR CONDUCTING AND SEMICONDUCTING SURFACES 


L INTRODUCTION 

With the attention toward increasingly lower energy (<50 eV) ion behavior in the 
magnetosphere,[ 1] there is a greater emphasis placed on the focusing properties of space plasma 
instruments designed to measure the flux of very low-energy ions. This need comes 
simultaneously as instrumentation development permits better resolution and with the capability 
to control spacecraft charging and differential charging. [2] Often in these instruments, the 
particle energies are selected electrostatically and their arrival angle, if determined, can be 
selected either electrostatically [3] or by collimation. [4] Using these techniques to measure 
low-energy particles, conditions arise that may cause particle trajectories to vary from their ideal 
path in the instrument. For example, the analysis of low-energy ions with retarding potential 
grids is often seriously affected by space charge and stray electric and magnetic fields. [5] 

The process that will be focused on in this paper is the effect of induced charges on 
low-energy particles as they travel near conducting and semiconducting surfaces. In particular, 
the question to be answered is the effect of the induced charge on these particles as they pass 
through grounded channels that are similar in dimension to those that may be designed in space 
plasma instruments. The Thermal Ion Dynamics Experiment (TIDE) on the future Global 
Geoscience (GGS) mission is designed to have an array of beryllium copper channels to serve as 
angle collimators for incident particles, each channel having a height of 1.27 mm, length of 25.4 
mm, and thus a length to width ratio of 20. Another example occurs in a design by Stenzel in 
which he replaces the first grid of a retarding potential analyzer with a grounded microchannel 
plate (MCP) resulting in a vary narrow angular response. The MCP is a honeycomb-like 
structure of hollow glass channels that are treated in such a way as to render the channel walls 
semiconducting with a resistivity ranging from 10 10 to 10 15 £2 cm. Typically these glass channels 
have internal diameters ranging from 0.01 to 0.1 mm and length to diameter ratios ranging from 
40 to 100. 

To answer the question posed above, a computer model was developed to inject singly 
charged particles into channels of both grounded, conducting plates, and semiconducting plates. 
The induced charge on the particles was then found by using the method of images. The 
dielectric constant for the semiconducting plates was e = 8.3, the same constant for the Coming 
8161 glass used in MCP plates provided by Galileo. [6] To study the effect of the induced 
charges various geometric parameters were varied: the parallel plates, ranged in length from 0. 1 
mm to 50 mm, the aspect ratio, defined here as the plate length to the separation of the plates, 
ranged from 1 to 100, and the distance of the simulated detector from the end of the plates varied. 



In the following, the design of the computer simulation will first be discussed, examples of 
particle trajectories with and without the force of induced charges will be shown, and results will 
be given on the effected channel geometric factor. 


0. PROCEDURE 


It is well known that for a point charge between parallel, grounded plates you can find the 
induced charge on the plates by the method of images. Although the problem can be solved by a 
wide variety of methods, [7-13] the method of images was chosen which is straightforward and 
can lead to an absolute and uniformly convergent series when treated carefully. 

The parallel planes were chosen to be at the locations of y = dfl andy = -d/2 and the 
charge was placed at an arbitrary distance, y, from the center of the channel (see Appendix A). 

To force the potential of the plane at y = d/2 to zero we place an image charge opposite in sign 
and symmetric to the plane, etc. Now to force the potential to zero on the plane at y = -d/2 we 
have to introduce two corresponding charges of opposite sign at the points symmetric to this 
plane, etc. An infinite series results and the series converges because each successive image 
decreases in magnitude. The resulting force on the charge q for conductive plates is 

lF , = iV y L , 

r-hy 0 2 - 4A 2 ) 

where j is tne image pair index, d is the distance between the plates, and h = y/d. The force on 
the charge q using semiconductive plates with a dielectric constant of 8.3 is 


lF/ = 



i 

(f-4h 2 ) 2 


(see Appendix A for the derivations). From these equations it can be seen that as j increases, the 
terms in the series decrease for a constant position of the particle. Also as seen from Appendix 
A, the even image pairs cancel so the sum is carried out over the odd values of j. For this study j 
was summed to 7 for the sake of computer efficiency. This introduced a percent error of less 
than 0.4 using an h of 0.01. Figure 1 shows a log-linear plot of the force as a function of h, 
showing the force on the particle is zero at the center of the channel, h = 0, and reaches infinity 
when the particle is near the plate at h = 0.5. This curve was derived using the force obtained 
from conductive plates. 

The computer model was constructed so that particles were injected between simulated 
parallel plates at incremental position and angle steps so as to illuminate the channel with all 
possible trajectories (see Figure 2). Transmission was then determined for the initial conditions. 
The angles were limited to the range that allowed trajectories with the greatest angle of 
incidence, this maximum angle defined by a = tan' 1 dJL t where L is the length of the channel. 
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The positions of entry, dy, were held between 0 and d/2 due to the symmetry of the problem. 
The maximum for both angles and positions was divided into steps of 24 creating a 24 x 24 dy- a 
transmission matrix. To illustrate, beginning at the first entry position at the center of the 
channel, a particle was injected at each of the ±12 angles, etc. The 25th particle then entered at 
the next entry position at the first angle, etc. The particles moved according to the equation of 
motion using the "leap frog" method in time steps resolved to 0.001 of their time to cross the 
channel. This resolution was a trade-off between accuracy of the trajectory and computer time. 
A simulated detector, the same width of the plate separation, was placed at the end of the 
channel. The distance of this detector from the end of the channel, or the force free region, was 
modeled in three different positions: flush to the channel or a distance of 1L or 2L from the end 
of the channel. Now with this simulation as just described, the parameters studied were the 
particle energy, the channel length, the aspect ratio, and the detector distance. 


ID. RESULTS 

Figures 3 and 4 illustrate the effect of induced charges on representative low-energy 
particles as they pass through grounded, conducting channels. In these figures, the vertical axis 
shows the width of the channel and the detector, and the horizontal axis shows the center of the 
channel and extends to the detector. The channels in these two figures have a width of 0.25 mm, 
a length of 2.5 mm, and the detector is placed one length of the channel away at 5 mm. Both 
panels of Figure 3 show the same five particles of 10 3 eV entering at the same positions near the 
plate and having the same incident angles of -5° to 6°. Panel a displays the trajectories with no 
induced force and panel b displays these trajectories with the induced force. Note that with the 
induced force, four out of the five particles are counted by the detector, whereas the particles 
with no force are not transmitted to the detector. Panels a and b of Figure 4 illustrate the induced 
force on particles of different energies. Again both panels show particles of the same entry 
positions, this time near the center of the channel, with the arrival angles -1° to 2.5°, but the 
particles in panel a are 1 eV and those of panel b are 10' 3 eV. As these particles pass through the 
channel you can see again, as in Figure 3, the stronger induced force as the particles near the 
plates. 


The nature of the trajectories, as affected by the induced charges, is more evident in 
Figures 5 and 6. Figure 5 shows dy - a transmission matrices whose elements are entry position 
and angle differentials with an * in the element marking the transmission of particles with 
corresponding conditions. The columns are labeled by entry element, dy, and the rows by angle 
element, a. These matrices were developed using once again a channel length of 2.5 mm, a 
channel aspect ratio of 10, and with the detector one length of the channel away at 5.0 mm. 
Panel a shows the transmission of 10‘ 2 eV particles with no induced force showing the 
characteristic window and panel b shows the transmission of 10' 2 eV particles with the induced 
force. Panel b illustrates that particles that entered near the plate had to enter at a more negative 
angle (away from the plate) to be transmitted. 
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Fig. 3. Examples of the effect of the induced force on low-energy particle trajectories within a 
channel. 
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Fig. 5. Transmission matrices. 
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Figure 6 shows phase plots developed from the transmission matrices similar to Figure 5. 
The outside dotted region shows the area created by the 24 x 24 position and angle elements. 

The single striped region is the area created by passing the particles through the channel with no 
induced force as seen in Figure 5a, and the region overlayed, shown with cross-hatches, is 
created by passing the same particles through the channel with the induced force as in Figure 5b. 
Both panels of Figure 6 were obtained with the same channel geometry as in Figure 5 but panel a 
has particle energies of 10 2 eV while panel b has particle energies of 10^ eV. Panel a shows a 
27% reduction in transmission while panel b shows a 95% reduction. 

To find the channel geometric factor as effected by the induced charge, we must know the 
relation of the number of particles CR (particles/accumulation) transmitted by the channel to the 
directional intensity, I(y,a,E), usually expressed in (particles/cm 2 sr eV accumulation). This 
relation is defined by 


CR = jl g(y, a, E)dAdQdE 

where y is the position vector, a is the angle of incidence, and E is the particle energy. g(y,a,E) is 
a function in the energy-dependent geometric factor which will be defined for one dimension as 

GF(E) = fg(y, a,E)dyda. 

This quantity is a function of the efficiency and geometry of the analyzer/detector and is 

expressed in units of cm 2 sr. For the study, the function, g, has the value of 1 for values of y and 
a corresponding to trajectories which strike the detector and zero otherwise. 

The 3D channel geometric factor can be obtained by resolving the particle motion due to 
the induced force into two perpendicular surfaces. Each calculation has the dimensions of length 
times angle, and then the total geometric factor is approximately the product of the two partial 
factors. This approximation is good for the small solid angles normally found in electrostatic 
analyzers. [14] 

Figure 7 shows the effect of the induced charges on the energy-dependent geometric 
factor. This geometry factor is normalized by the geometric factor obtained using no induced 
force. These values were calculated using conductive plates. The dimensions of the channels 
used in both 7a and 7b have an aspect ratio of 10 but differ in that panel a has a channel length of 
2.5 mm, while panel b has a channel length of 0.25 mm. The three curves from right to left mark 
the distance of the detector from the channel aperture; flush, 1 L, and 2 L, respectively. It can be 
seen by comparing both panels a and b that as the channel length decreases by a factor of 10 the 
particle energy has to increase by a factor of 10 to have the same geometric factor. Also note 
that at 1 eV the geometric factor is already reduced by 10%. 

As a measure of the perceptibility of the induced charge effect, the value of the particle 
energy where the geometric factor fell by 50% was then plotted versus channel length ranging 
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Fig. 7. Normalized geometric factor as a function of particle energy. 
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from 0.1 to 50 mm in Figure 8. This was done for both conducting and semiconducting channels 
with aspect ratios of 10. The energy value, denoted by ^0.5’ is proportional to the length of the 
channel, L, by a constant and is defined as: 


E 0.5 


C 

T 


where C = 0.006 eV mm for conductors and 


Eq.5 


Cd 

= X 


where C D = 0.0046 eV mm for a semiconductor with a dielectric constant of 8.3. The detector is 
placed at 1 L for each channel length where L is in mm. 

Figure 9 shows the energy values plotted as a function of aspect ratio for both conductors 
and semiconductors. A relation between ^0.5 and channel dimension was found and can be used 
for any arbitrary channel. This relation is defined as 

£ C ( aspect ratio ^ 

where C = 0.006 eV mm for conductors and 

„ Cd ( aspect ratio \ 3 

Eo.s =77 V r iff ) 


where C D = 0.0046 eV mm for semiconductors. 


IV. CONCLUSIONS 

The induced force was found on particles as they passed through grounded, parallel plates 
which served as channels ranging in length from 0.1 to 50 mm and in aspect ratio from 1 to 100. 
This force was calculated for both conducting and semiconducting plates. The reduced 
geometric factor was found as a function of particle energy as the channel length, aspect ratio, 
and distance of the detector from the end of the channel was varied. Relations were obtained 
between the channel dimensions and the energy value where the geometry factor fell by 50%. 

For comparison, the angle collimators for TIDE have lengths of 25.4 mm and aspect ratios of 20. 
From the equations above we know that for conducting channels that E os will be 0.004 eV. For a 
grounded MCP with channel diameters of 50 pm and a length to diameter ratio of 60, E Q 5 is 0.3 
eV. The results from this study show a need to consider the effect of induced charges on 
low-energy particles less than 1 eV as they near conducting or semiconducting surfaces for the 
channel dimensions cited. Although the channels of MCPs are weakly conducting, for these low 
energies and their small channel dimensions the induced force on the particle does need attention. 
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As a rule of thumb, the effect is appreciable at the 0.1 eV level for aspect ratios larger than 25 
and lengths shorter than 1 mm. 
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Appendix A. Derivation of the Induced Force on a Charge Between Two Plates 
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The plates coincide with the xz plane a distance d apart, and the point charge is at an arbitrary 
distance y above the center axis. The distance between the point charge and each image charge is 
r t] where i is the number shown on the plates and j is the image pair index. The distance from the 

point charge of the first few images and their polarities are given below: 
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Note that every even image pair cancels. 
Using the equation 


F = 


_ _v 


4k£ 0 


oo __ _ 

Y J L 

l_ r tj 2 r 2 i*J 


the resulting force on the charge, < 7 , using conductive plates is 


F = 


_ V. 


4jt£ 0 


jjfes... _(oV/-4y) 2 ) ' ((/d+4y) 2 )] 


where j is the image pair index. Simplifying, this equation becomes 


Ml 


F _ +£_ y 

" Ms... if* -*ft 


F = 


+2 <?h y j 

Ms... if - 4 A 2 ) 2 


y 

where h = jj. 


To calculate the force on the charge, q, using dielectric boundaries, the image charges are 
replaced by q ('e ; -e 2 /e ; + e 2 j[15] where e, is the permittivity within the channel and ^ is the 

permittivity of the dielectric boundary. The resulting force is 

r <f ( Ei -£2 ^ y %jdy 
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